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ABSTRACT 

We have detected the Rossiter-Mclaughlin effect during a transit of WASP-47b, the only known hot Jupiter 
with close planetary companions. By combining our spectroscopic observations with Kepler photometry, we 
show that the projected stellar obliquity is A = 0° ±24°. We can firmly exclude a retrograde orbit for WASP- 
47b, and rule out strongly misaligned prograde orbits. Low obliquities have also been found for most of the 
other compact multiplanet systems that have been investigated. The Kepler-56 system, with two close-in gas 
giants transiting their subgiant host star with an obliquity of at least 45°, remains the only clear counterexample. 


Subject headings: planetary systems — planets and 
rotation 

1. INTRODUCTION 

Stars that host close-in giant planets have been found to 
have a wide range of obliquities (see, e.g., Triaud et al. 2010, 
Albrecht et al. 2012, or the review by Winn & Fabrycky 
2015). Attempts to explain this surprising result have until 
recently focused on the implications for the formation mech¬ 
anism of hot Jupiters. The broad range of obliquities seemed 
to suggest that whichever processes shrink the orbits of hot 
Jupiters are also capable of tilting their orbits away from the 
plane of formation. It remains possible, though, that the obliq¬ 
uity distribution has nothing to do with hot Jupiters per se, but 
is instead the outcome of more general post-formation dynam¬ 
ical processes. 

This ambiguity arises because measurements of stellar 
obliquity have been largely confined to hot Jupiters, for prac¬ 
tical reasons. The main technique used to measure obliquities 
is the Rossiter-McLaughlin (RM) effect, which is easier to ob¬ 
serve for larger planets in short-period orbits. More recently it 
has become possible to probe the obliquities of other types of 
systems, particularly the compact systems of multiple transit¬ 
ing planets found using data from the Kepler spacecraft. The 
RM effect is more difficult to observe in such systems because 
the planets are generally smaller and the transits are less fre¬ 
quent than for hot Jupiters, but detections have been achieved 
for two such systems; Kepler-25 and Kepler-89 (Hirano et al. 
2012, Albrecht et al. 2013). Both stars have low obliquities. 
Low obliquities were also measured for the Kepler-30 sys- 
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tern by analyzing starspot-crossing events (Sanchis-Ojeda et 
al. 2012) and for the Kepler-50 and Kepler-65 systems by an¬ 
alyzing the rotational splittings of p-mode oscillations (Chap¬ 
lin et al. 2013). Thus it seemed that the stellar hosts of com¬ 
pact multiplanet systems generally have low obliquities (Al¬ 
brecht et al. 2013X and that the high obliquities were indeed 
peculiar to hot Jupiters. 

This tentative conclusion was shattered by the discovery 
that the host star of Kepler-56, a compact system of two tran¬ 
siting planets, has an obliquity of at least 45° (Huber et al. 
2013). Batygin (2012) anticipated the existence of such sys¬ 
tems, as the outcome of disk migration in the presence of a 
distant companion star on a misaligned orbit. Li et al. (2015) 
demonstrated that the misalignment of Kepler-56 could have 
also arisen independently of disk migration, as a consequence 
of planet-planet dynamics. Thus, it now seems possible that 
the obliquities of hot-Jupiter hosts gave us a glimpse of more 
widespread phenomena, and that we might expect a wide 
range of stellar obliquities for planetary systems in general. 

In this light, the WASP-47 system has a special importance, 
as the only known case of a hot Jupiter that is also part of 
a compact multiplanet system. The hot Jupiter in this “hy¬ 
brid” system has an orbital period of P = 4.1 days and was 
discovered several years ago (Hellier et al. 2012). Subse¬ 
quently, long-term radial-velocity monitoring revealed a Jo¬ 
vian planet with a much longer period (Neveau-VanMalle et 
al. 2015, submitted). Most recently, data from the K2 mis¬ 
sion (Howell et al. 2014) revealed that there are two addi¬ 
tional transiting planets, a super-Earth with a period of 0.79 
days and a Neptune-size planet in a 9-day orbit (Becker et al. 
2015). Such configurations are thought to be rare (Steffen et 
al. 2012). Measuring the obliquities of the stars in such sys¬ 
tems might help to determine the relationship between them, 
the other compact multiple-planet systems, and the more iso¬ 
lated hot Jupiters. 

This paper is organized as follows. Section |2] presents the 
spectroscopic and photometric observations of WASP-47 that 
were used in our study. Section[3]presents our analysis lead¬ 
ing to the determination of the projected stellar obliquity. Sec- 
tionlUdiscusses these results and compares WASP-47 to other 
known planetary systems. 
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2. OBSERVATIONS AND DATA REDUCTION 

WASP-47 happened to fall within the field of view of 
“Campaign 3” of the K2 mission, which uses the Kepler space 
telescope to perform precise photometry of stars in selected 
ecliptic fields (Howell et al. 2014). Becker et al. (2015) have 
presented an analysis of the photometric data, including the 
discovery of the additional planets “d” and “e”. Here we use 
the short-cadence photometric time series provided by Becker 
et al. (2015), which have a time sampling of one minute and 
were corrected for long-term trends due to stellar variability 
and instrumental artifacts. 

In an attempt to detect the RM effect, we monitored the 
optical spectrum at high resolution using the HIRES spectro¬ 
graph (Vogt et al. 1994) on the Keck I telescope, on the night 
of UT August 10, 2015. The weather was excellent, with sub- 
arcsecond seeing. We obtained 29 spectra over the course of 
5 hours spanning a transit of the hot Jupiter WASP-47b. Ex¬ 
posure times were 11-13 minutes, set by an exposure meter 
that enforced a signal-to-noise (S/N) ratio of 100 pixeE^ near 
550 nm for each spectrum. We used the 0.87"-wide C2 decker 
giving a spectral resolution of 55,000. A glass cell of molecu¬ 
lar iodine was placed in the light path to calibrate the HIRES 
spectra. The next night, we obtained a spectrum of WASP-47 
without the iodine cell, and with S/N = 200 pixeP*, to serve 
as the template spectrum in the Doppler analysis. 

The reduction of HIRES spectra and the Doppler analy¬ 
sis used the standard pipeline of the California Planet Sur¬ 
vey (Howard et al. 2010). The deconvolved iodine-free spec¬ 
trum of WASP-47 was used with the molecular iodine line 
atlas in a forward model that determines the wavelength zero 
point, spectral dispersion, Doppler shift, and the instrumen¬ 
tal profile, following the methods of Marcy & Butler (1992). 
The forward model works on 2 A spectral segments and uses 
the weighted average to determine the final radial velocity for 
each of the observations. The results are shown in Figure [T] 
and presented in Table[T] 

3. ANALYSIS OE THE OBSERVATIONS 
3.1. Determination of projected rotation rate 

The transits of WASP-47b are very nearly central, i.e., the 
transit impact parameter b is close to zero. In such situa¬ 
tions, analyses of the RM effect are subject to a strong de¬ 
generacy between the sky-projected stellar obliquity A and the 
sky-projected stellar rotation rate v sin i (Gaudi & Winn 2007). 
Therefore it is of special importance to determine or constrain 
the value of v'sin/ independently of the RM effect. 

We analyzed the iodine-free, high-S/N spectrum using 
SpecMatch (Petigura 2015). This is a code for obtaining the 
basic spectroscopic parameters from high-resolution spectra. 
It creates synthetic spectra for any input values of Teff, log.g, 
and [Fe/H], by interpolating between a grid of model spectra 
taken from the library of Coelho et al. (2005). SpecMatch 
also applies broadening kernels to account for rotational and 
macro-turbulent broadening. Importantly, SpecMatch also 
measures the HIRES instrumental line broadening, which de¬ 
pends on the slit width and atmospheric seeing at the time of 
observations, by fitting the O 2 B-X telluric lines with a comb 
of Gaussians. For WASP-47 the instrumental profile is mod¬ 
eled as a Gaussian function with a dispersion of 1.95 km s“'. 

Given an observed spectrum, SpecMatch adjusts the 
spectroscopic parameters of the synthetic spectrum until it ob¬ 
tains the best match. For WASP-47, the results were Teff = 
5565 ± 60 K, logg = 4.29 ± 0.07, and [Fe/H] = 0.43 ± 0.04. 


TABLE 1 

Radial velocity observations 


Time [BID] 

Radial velocity tm/s] 

Uncertainty tm/s] 

2457244.871873 

24.1 

1.4 

2457244.879755 

29.7 

1.6 

2457244.887822 

22.8 

1.5 

2457244.896063 

19.5 

1.7 

2457244.904292 

20.7 

1.6 

2457244.912684 

20.0 

1.6 

2457244.920959 

17.6 

1.5 

2457244.929316 

16.2 

1.7 

2457244.937406 

14.6 

1.6 

2457244.945624 

20.8 

1.6 

2457244.953888 

22.4 

1.5 

2457244.962233 

22.6 

1.6 

2457244.970729 

18.9 

1.6 

2457244.979537 

12.9 

1.6 

2457244.988021 

4.2 

1.6 

2457244.996621 

0.7 

1.6 

2457245.005799 

-6.0 

1.7 

2457245.014526 

-9.0 

1.6 

2457245.023033 

-16.5 

1.5 

2457245.031413 

-17.7 

1.6 

2457245.039723 

-20.2 

1.5 

2457245.048254 

-25.4 

1.6 

2457245.057004 

-29.9 

1.5 

2457245.066506 

-29.9 

1.8 

2457245.075314 

-24.3 

1.5 

2457245.083902 

-21.3 

1.5 

2457245.092618 

-24.9 

1.5 

2457245.101518 

-33.1 

1.5 

2457245.110824 

-29.9 

1.5 


These values are all in agreement with those presented by 
Hellier et al. (2012) and Mortier et al. (2013). For the 
projected rotation rate, SpecMatch found an upper bound 
V sin i = 0 ± 2 km s“', based on the lack of detectable rotational 
broadening. This result is slightly different from the value of 
3.0 ±0.6 km s“' reported by Hellier et al. (2012). To check 
on the SpecMatch result, we also performed an independent 
analysis of the line profile using the tools described by Hirano 
et al. (2014). We found vsin/ = 1.3];[ 2 km s“', in good agree¬ 
ment with the value found by SpecMatch. In the subsequent 
analysis we use vsin; = 0±2 km s“', because it is based on 
a higher-S/N and higher-resolution spectrum, and because we 
have calibrated the vsin / scale of SpecMatch against all our 
previous RM results. We checked that this choice did not af¬ 
fect our final conclusions significantly, by using the other two 
constraints, obtaining similar results. 

3.2. Determination of the projected obliquity 

Our first goal was to create a phase-folded light curve for 
the transits of WASP-47b, which could be fitted simultane¬ 
ously with the Doppler data. For this purpose we used the 
detrended photometric time series provided by Becker et al. 
(2015). We generated an initial phase-folded transit light 
curve for WASP-47b, using the ephemerides of Becker et al. 
(2015) to fold the data with the constant period of WASP-47b 
and discarding the data that were obtained during transits of 
planets d and e. 

Next, we wanted to improve on this light curve by tak¬ 
ing into account the transit-timing variations of WASP-47b. 
We fitted the initial light curve with the model of Mandel 
& Agol (2002) for the case of a quadratic limb-darkening 
law. The free parameters were the transit depth 
impact parameter b, the scaled stellar radius and the 

limb-darkening coefficients ui and M 2 - Then the best-fitting 
parameters were held fixed, and the model was used as a tern- 
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Fig. 1.— RM model fit to the radial velocities during a transit of 
WASP-47b. Radial velocity observations are represented with black dots, 
with uncertainty bars that are sometimes smaller than the dots (and do not 
include the jitter term). The best-fit RM model to the data shown in the upper 
panel gives a A near zero, and it is represented with a red thick line. The 
lower panel shows the residual radial velocity respect to the model. 

plate to determine the midtransit time of each individual tran¬ 
sit. We fitted the 5 hours of data spanning each transit with 
3 free parameters; the midtransit time, and the two parame¬ 
ters describing a linear function of time (to take into account 
the covariance between the midtransit time and longer-term 
trends in stellar brightness). The series of midtransit times de¬ 
termined in this manner were well-fitted by a quadratic func¬ 
tion of transit epoch, which takes into account the observed 
transit timing variations that have a maximum deviation of 1 
minu tf3- This new orbital ephemeris was then used to re¬ 
fold the original time series, producing a new phase-folded 
light curve for WASP-47b that takes into account the small 
transit-timing variations. To speed subsequent computations 
we averaged the light curve into one-minute bins. The result¬ 
ing light curve is well-fitted by a standard transit model, with 
residuals that appear random and independent. To determine 
appropriate weights for the binned data points, we assigned 
an uncertainty to each point such that x^in ~ -^dof- 

We then fitted simultaneously the radial-velocity data and 
the binned photometric data. To describe the photometric 
transit we used the same Mandel & Agol (2002) model with 
the same five free parameters, but with both limb-darkening 
coefficients subject to Gaussian priors based on theoretical es¬ 
timates (Claret & Bloemen 2011). The mid-transit time of the 
folded light curve was allowed to vary freely. 

To describe the RM effect, we used the code of Hirano et 
al. (2011). This code takes into account all the important line¬ 
broadening mechanisms as well as the RM effect. The key 
parameters were vsini and A. For vsin/, we required it to 
be positive, and used a Gaussian prior with a mean of zero 
and a dispersion of 2 km s“*, based on our findings from the 
previous section. We also modeled the out-of-transit radial- 
velocity variation with a linear function of time, giving two 
additional parameters 7 and 7. We adopted a Lorentzian (nat¬ 
ural) broadening of 1 km s“* and a Gaussian (instrumental 
plus thermal) broadening of 3.3 km s“', based on the instru¬ 
mental profile of HIRES and the effective temperature of the 
star. The macro-turbulent broadening was allowed to vary, 

*** Becker et al. (2015) modeled the timing variations as a sinusoidal func¬ 
tion of epoch, which gives essentially the same results. 


TABLE 2 

Key parameters of WASP-47b 


Parameter 

Value 

68.3% Conf. Limits 

Sky-project obliquity, A [deg] 

0 

±24 

vsin / [km s“*] 

1.80 

(1-0.24,-0.16) 

RV offset, 7 [m s“*] 

-4.3 

±0.7 

RV slope, 7 [m s“* day“*] 

-244 

±8 

Stellar jitter, a [m s“^ ] 

2.1 

±0.5 

Midtransit time [BJDtdb] 

2457245.0049 

±0.0017 

iRp/R*f 

0.01034 

±0.00004 

Ri,/a 

0.1025 

±0.0007 

Impact parameter, b 

0.13 

(1-0.05,-0.08) 


subject to Gaussian prior corresponding to 3 ±0.5 km s“*, 
based on the relationships given by Doyle et al. (2014). The 
convective blueshift was taken into account with an assumed 
amplitude of 0.5 km s“', following Albrecht et al. (2012). We 
decided to allow the midtransit time to vary freely, rather than 
constraining it to conform to the transit ephemeris, given the 
observed level of transit-timing variations in the system. 

The last parameter of the model is the “stellar jitter” term, 
which is added in quadrature to the radial-velocity uncertain¬ 
ties that are internally estimated by our Doppler code. We 
decide to allow the jitter term to vary freely, so rather than 
minimizing the function, we minimized the negative of the 
logarithm of the likelihood. In this way we can naturally con¬ 
trol the jitter term by adding the log of the final radial velocity 
uncertainty (see equation 1 of Johnson et al. 2011). This en¬ 
sures a similar result to what we would find fixing the jitter 
term requiring the best-fit x^ = Ndof, 

We determined the allowed ranges for the 13 model pa¬ 
rameters, and their covariances, with a Monte Carlo Markov 
Chain routine. Figure [T] shows the best-fitting model to the 
RV signal. The results for the marginalized distribution of 
each relevant parameter can be found in Table |2] Our pho¬ 
tometric parameters are all compatible with those found by 
Becker et al. (2015), including the midtransit time. The stel¬ 
lar jitter of 2.1 m s“' is a typical value for a star of this type. 
From the value of 7 we can get an estimate of the radial ve¬ 
locity semi-amplitude K of ITrm =-jPorb/'2.Tr = 161 ±5 m s“*, 
slightly larger than the value ofIir=136±5ms“' found in the 
discovery paper (Hellier et al. 2012). 

As expected the uncertainty in A is dominated by the co- 
variance with vsin/, a characteristic of systems with a low im¬ 
pact parameter. Figure|2]shows the joint posterior distribution 
function of A and vsin / resulting from our MCMC analysis. 
The covariance between A and vsin / causes the uncertainty in 
A to be larger than one might expect from such a high-S/N de¬ 
tection of the RM effect. If it were possible to determine v sin / 
accurately and precisely from external observations, this de¬ 
generacy could be broken, and a precise value of A could be 
obtained, but this is not yet the case for WASP-47. 

4. DISCUSSION 

We have presented the third measurement of the RM effect 
for a compact multiplanet system. In all three cases the pro¬ 
jected obliquity has been found to be consistent with zero. For 
WASP-47 our determination of A was limited by the strong 
covariance between A and vsin/, an effect caused by the low 
impact parameter of the planetary transits. In principle this 
problem could be ameliorated by a better determination of 
the impact parameter, but this is unlikely to be achieved soon, 
given that the K2 data already has such high photometric pre- 
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Fig. 2.— The degeneracy between A and vsin; increases the uncer¬ 
tainty on A. The central plot shows a contour plot for the combined distribu¬ 
tion of values for A and v sin i from our MCMC routine. The three solid lines 
represent the 1, 2 and 3cr contours. Darker regions represent areas with higher 
posterior density. A better knowledge of v sin i (beyond our upper bound of 
2 km s“') would allow a better determination of A. 


cision. Another way to improve the measurement of A would 
be to refine our constraint on vsin/, either via new and higher- 
resolution spectroscopic observations, or by measuring the ro¬ 
tation period of the host star. Regarding the latter, Hellier et 
al. (2012) did not report any quasiperiodic variability in the 
WASP data that could be attributed to stellar rotation, and 
likewise, we have not found any convincing signals in the K2 
data. 

Our results have allowed us to exclude a retrograde orbit 
for WASP-47b, and to rule out strongly misaligned prograde 
orbits. As recently as 2008, this would not have been the 
least bit surprising. In particular, the theory of disk migra¬ 
tion — a leading contender to explain the existence of close- 
in planets and multiplanet systems — was expected to lead 
to well-aligned and low-obliquity systems. But in the past 
few years we have learned that stellar obliquities can vary 
widely, for reasons that are not yet understood. For WASP- 
47, a low obliquity or even a moderate obliquity was not a 
foregone conclusion, especially given how unusual the sys¬ 
tem is in comparison to the other known systems. 

Consider, for example, a comparison between WASP-47 
and other systems discovered by Kepler with at least one tran¬ 
siting planet with a period shorter than one day (the “ultra- 
short-period” planets studied by Sanchis-Ojeda et al. 2014). 
The Kepler telescope has revealed additional transiting com¬ 
panions in these systems, but these companions are typically 
Neptune-sized or smaller. In only one case, KOI 191, is there 
a Jupiter-sized companion planet similar to WASP-47b, and 
in that system the giant planet has a period of 15 days. 

Or consider the class of systems with one close-in giant 
planet and at least one additional planet. Using the online 
planet catalog exoplanets.org (Han et al. 2014), we 
found that among the 205 planetary systems with a transit¬ 
ing planet larger than 0.5 Rjup and orbital periods shorter 
than 15 days, only 5 of them have known planetary compan¬ 
ions. Some of these are Jupiter-size planets with very long- 
period companions (see, e.g., Knutson et al. 2014). Others 
have planets that are half the size of Jupiter and have smaller, 
shorter-period companions (such as Kepler-18, Cochran et al. 


2011). None of them resemble WASP-47, which has a radius 
of 1.1 Rjup, as well as nearby companions both interior and 
exterior. 

Finally, consider a comparison with a mass-selected sam¬ 
ple. There are 245 known stars with at least one planet with 
a mass larger than 0.1 Mjup and a period shorter than 15 days. 
Only 11 of these have more than one known planet, many 
of which were also in the radius-selected sample described 
above. The only system within this sample that resembles 
WASP-47 in having a close-in giant planet with a small inte¬ 
rior companion and at least one exterior companion is 55 Cnc 
(Butler et al. 1997, Marcy et al. 2002, Fischer et al. 2008, 
Winn et al. 2011). This planetary system has a giant planet 
with an orbital period of 14.6 days, and a 2 Rg planet with a 
period of 17.7 hours, along with a few other gas giants with 
periods of 44, 260 and 5170 days. Thus, WASP-47 might be 
considered a more compact cousin of the 55 Cnc system. In 
this respect it is interesting to note that Bourrier & Hebrard 
(2014) found 55 Cnc to have a projected obliquity of 72°, al¬ 
though this conclusion has been disputed by Lopez-Morales 
et al. (2014). 

We would need to extend the orbital period range of the 
transiting Jupiter-size planet up to 25 days to find another ex¬ 
ample of a system with a planet smaller than Neptune interior 
to the orbit of a warm Jupiter-size planet. The Kepler-89 sys¬ 
tem (Weiss et al. 2013) contains a mini-Neptune planet with 
an orbital period of 10 days, a gas giant with a period of 22.3 
days and an inner smaller planet with an orbital period of 3.7 
days. In this case there are two planets interior to the orbit of 
the Jupiter-size planet, but they are not as close to the host star 
as the super-earths 55 Cnc e and WASP-47d. 

It remains the case that the only clear and undisputed exam¬ 
ple of a compact multiplanet system with a high stellar obliq¬ 
uity is Kepler-56 (Huber et al. 2013). Additional measure¬ 
ments of obliquities in multiplanet systems, either in unusual 
systems such as WASP-47 or using other techniques more ca¬ 
pable to study larger samples of diverse planet hosts (Morton 
& Winn 2014; Mazeh et al. 2015), may play an important 
role in future attempts to understand their formation and their 
relationship to more isolated hot-Jupiter systems. 
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